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FIELD   DETERMINATIONS   OF    A   PROBABILISTIC   DENSITY    FUNCTION 
FOR    SLOPE    STABILITY   ANALYSIS   OF   TAILINGS    EMBANKMENTS 

By  P.  C.  McWilliams1   and  D.  R.  Tesarik2 


ABSTRACT 

Theoreticians  in  soil  mechanics  have  been  pursuing  a  probabilistic 
approach  to  the  factor  of  safety  of  an  embankment  or  dam  for  the  past 
10  years.  The  motivation  for  this  work  is  in  contrast  to  the  current 
practice  of  deterministically  computing  a  factor  of  safety  and  treating 
it  as  an  absolute  with  no  regard  to  its  inherent  statistical 
variability. 

Basic  to  the  probabilistic  approach  is  the  selection  of  an  appropri- 
ate statistical  model  to  represent  the  histogram  or  probability  density 
function  (PDF)  of  the  factor  of  safety  values.  Rather  than  simply  as- 
suming which  PDF  is  appropriate,  the  Bureau  of  Mines  collected  data  at 
two  waste  disposal  embankments  for  consideration.  This  report  ad- 
dresses three  candidate  models,  using  the  techniques  of  nonlinear  curve 
fitting,  and  identifies  the  "best"  model.  A  propagation  of  error  for- 
mula for  estimating  the  variability  of  Fellenius'  factor  of  safety  is 
also  discussed. 


Mathematical  statistician. 


1 

2Mathematician. 
Both  authors  are  with  the  Spokane  Research  Center,  Bureau  of  Mines,  Spokane,  Wash. 


INTRODUCTION 


Theoreticians  in  soil  mechanics  have 
been  pursuing  a  probabilistic  approach 
to  the  factor  of  safety  of  an  embank- 
ment or  dam  for  the  past  10  years.  The 
motivation  for  this  work  is  in  contrast 
to  the  current  practice  of  deterministi- 
cally  computing  a  factor  of  safety  and 
treating  it  as  an  absolute  with  no  re- 
gard to  its  inherent  statistical  vari- 
ability. Stated  simply,  in  the  past, 
soils  engineers  have  used  a  single  pa- 
rameter, the  factor  of  safety.  It  would 
now  seem  prudent  to  impose  an  additional 
parameter,  the  dispersion  of  the  fac- 
tor of  safety.  Given  these  two  param- 
eters, one  may  then  compute  the  proba- 
bility of  failure  as  a  measure  of  slope 
stability. 


Basic  to  the  probabilistic  approach  is 
the  selection  of  an  appropriate  statisti- 
cal model  to  represent  the  histogram  or 
probability  density  function  (PDF)  of  the 
factor  of  safety  values.  It  is  certainly 
easy  to  hypothesize  such  a  model;  the 
normal  PDF  is  almost  always  the  first 
choice.  For  this  reason,  the  Bureau  of 
Mines  collected  data  at  two  waste  dis- 
posal embankments  such  that  their  PDF's 
could  be  investigated.  This  report  is 
addressed  to  investigating  candidate  mod- 
els via  the  techniques  of  nonlinear  curve 
fitting,  and  then  choosing  the  "best" 
model  available.  A  propagation  of  error 
formula  for  estimating  the  variability 
of  Fellenius'  factor  of  safety  is  also 
discussed. 
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testing  of  this  large  set  of 


PROBABILISTIC  APPROACH 


There  are  several  prior  discussions  of 
applying  probability  theory  to  slope  sta- 
bility analysis  (l_-2,  ]_,  .12,  15-16, ).3 
For  orientation  purposes,  it  would  seem 
prudent  to  briefly  sketch  the  problem. 
Most  dam  design  engineers  use  the  methods 
of  either  Bishop-Morgenstern  or  Fellenius 
(_5_,  7_,  9)  to  compute  the  factor  of  safety 
for  slope  stability  design.  These  tech- 
niques are  deterministic  in  that  one  val- 
ue is  computed  which  measures  the  "safety 
of  the  dam."  The  basic  postulate  is  that 
of  a  potential  slip-circle  failure.  In 
terms  of  moments,  the  factor  of  safety  is 
the  ratio  of  the  resisting  moment  to  the 
overturning  moment  of  the  embankment 
(fig.  1).  If  this  ratio  is  greater  than 
1,  the  structure  is  assumed  to  be  in  a 
nonfailure  state.  Often  there  are  so- 
called  rules  of  practice;  for  example, 
a  factor  of  safety  of  1.5  or  greater  is 
required  in  some  applications. 

3Underlined  numbers  in parentheses  re- 
fer  to  items  in  the  list  of  references 
preceding  the  appendix. 


The  probabilistic  approach  suggests 
that  in  addition  to  computing  the  fac- 
tor of  safety,  one  should  also  include 
its  variability.  One  such  approach  is 
to  calculate  a  confidence  interval  for 
the  factor  of  safety  by  collecting  sev- 
eral samples  and  applying  the  theorem 
that  a  mean  value  (in  this  case,  the 
average  factor  of  safety  value)  is 
normally  distributed  (8).  This  is  dis- 
cussed in  detail  in  the  section  entitled 
"Confidence  Interval  for  the  Factor  of 
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Safety."  A  necessary  ingredient  in  most 
probabilistic  work  is  the  assumption  of  a 
PDF  that  characterizes  a  population  of 
factor  of  safety  values. 

Another  approach  is  to  compute  the 
probability  of  failure,  which  is  repre- 
sented by  the  area  under  the  PDF  to  the 
left  of  the  value  1  (fig.  2).  The  main 
objective  of  this  report  was  to  collect 
enough   sample  values   from  two  tailings 


embankments  to  establish  the  PDF  for  the 
factor  of  safety  values.  In  practice,  it 
would  not  be  expedient  to  expend  the  time 
and  money  to  collect  the  necessary  number 
of  samples  required  to  derive  the  PDF. 
Work  such  as  that  done  herein  should  pro- 
vide some  assurance  that  it  would  not  be 
unreasonable  to  postulate  a  particular 
PDF  without  actually  collecting  a  partic- 
ularly large  sample  of  data. 


PROPAGATION  OF  ERROR4 


An  interesting  alternative  to  collect- 
ing enough  samples  to  establish  a  PDF 
is  to  apply  the  statistical  technique 
known  as  propagation  of  error  (6).  This 
technique  will  approximate  the  standard 
deviation  of  the  factor  of  safety  (sp). 
It  may  not  be  convenient  for  the  designer 
to  collect  enough  data  to  compute  Sp  di- 
rectly; he  or  she  may  instead  assume 
values  for  the  coefficient  of  variation 
for  the  soil  properties — and  then  apply 
propagation  of  error.  Reference  is 
invited  to  either  Lee  and  Singh  (11)  or 
table  1  of  this  paper  for  candidate  co- 
efficient of  variation  values. 


Mathematical 
follows: 


considerations 


are  as 


Given  a  function  of  several  variables, 


G  =  G  (a,  b,  c), 


(1) 


4There 
between 


exist 
the 


notational 
disciplines 


ambiguities 
of    the 


mathematicians-statisticians  and  that  of 


the  geotechnical  engineers, 
per,  deference   is  given  to 


In  this  pa- 
the  notation 


of  mathematics  and  statistics. 


KEY 

F  =  Factor  of  safety 

F  =  Mean  factor 

of  safety 
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FIGURE  2.  -  Concept  of  probability  of  failure, 


by  propagation  of  error,  the  variance  of  G  (sG),  is  given  by 


Sg  =  Ga  sa  +  Gb  sb  +  Gc  sc  +  2rab  Ga  Gb  sa  sb 
+  2rac  Ga  Gc  sa  sc  +  2rbc  Gb  Gc  sc  sb, 
where  G\      =   partial  derivative  of  G  with  respect  to  the  i"^  parameter, 


(2) 


and 


S|   =  variance  of  the  i—  parameter, 

r|j  =  correlation  coefficient  for  i—  and  j—  parameters 


Assume  that  the  embankment   of  consideration  is  not  "zoned";   that  is,   it  is   of  a 
homogeneous  nature. 


Fellenius'  equation  for  the  factor  of  safety  (F)  for  soils  may  be  written  as 

-  cL  +  tan  <j)  £  (y  A|  cos  0[  -  U]  A£|) 
Y  Z|Aj  sin  0| 

where   c  =  cohesion  of  soil, 

<j)  =  angle  of  internal  friction, 

L  =  total  arc  length  of  slip  circle, 

Y  =  density  of  soil, 

Uj  =  pore  pressure,  i—  slice, 

Aj  =  area,  i—  slice, 

0j  =  angle  of  tangent  to  arc,  it-^  slice, 

AA|  =  arc  length,  i—  slice, 

and     £  =  summation  over  the  "i"  slices,  i  =  1,  n. 

TABLE  1.  -  Basic  statistics  for  soil  parameters 


(3) 


Slope  stability 
parameter 


Mean 


Cv  ,  pet 


Slope  stability 
parameter 


Mean 


Cv ,  pet 


CDA  SITE  (74  SAMPLES) 


SW  SITE  (130  SAMPLES) 


* ■ 

c 

Y 

F-Bishop..., 
F-Fellenius, 


-deg 
.psi 
.pcf 


36.74 

3.52 

114.00 

1.54 

1.43 


2.18 

2.94 

6.06 

.41 

.40 


5.9 
83.5 

5.3 
26.6 
28.0 


* 

c 

Y 

F-Bishop. .. 
F-Fellenius , 


.deg 
.psi 
.pcf 


39.13 

4.32 

106.85 

2.07 

1.95 


4.14 

4.14 

5.65 

.27 

.24 


10.6 
95.8 
5.3 
13.0 
12.4 


<}>  =  angle  of  internal  friction. 

c  =  cohesion. 

Y  =  soil  density. 

F  =  factor  of  safety. 


Mean  =  sample  average  value. 

S    =  sample  standard  deviation. 

Cw   =  coefficient  of  variation  = 


100-S 

■ 

mean 


See  figure  1  for  the  geometry  of  Fellenius'  method. 

Applying  the  propagation  of  error  formula  to  equation  3  yields 


,2   2 


sF  =  *>  s^  +  Fc  sc  +  Fy  sy  +  I   FU|  su 


2rc  ,yFcFYscsY  + 


.   +  2r 


un-1 


+   2rc,<|>FcF<|>scs<j> 


"n-run&un-run. 


(4) 


The  required  partial  derivatives  are 


u? 


and 


Y  E  A|  sin  G| ' 

sec2<{>  Z  (y  A|  cos  0|-  U|  &l\) 
~~y      £  A|  sin  0 

-tan  <}>  A£j 

Y  E  A|  sin  Oj ' 

-cL  +  tan  <f»  Zu|  Mj 


Y2  I  A|  sin  0, 


(5) 
(6) 
(7) 
(8) 


Of  particular  interest  is  the  selection 
of  Fellenius'  equation  rather  than  that 
of  Bishop.  Fellenius1  equation  gives  a 
"closed  solution"  for  the  factor  of  safe- 
ty value,  while  Bishop's  equation  is  not 
nearly  so  convenient,  with  an  iterative 
process  being  required  for  solution.  Al- 
so, many  practioners  and  soil  mechanics 
text  authors  still  seem  to  prefer  Felle- 
nius' formulation.  Later  in  this  report, 
the  propagation  of  error  methodology  will 
be  applied  to  the  collected  data  sets. 

A  proper  measure  of  the  uncertainty  of 
the  factor  of  safety  involves  more  than 


simply  computing  its  standard  deviation. 
Pore  pressure  uncertainty  is  a  very  im- 
portant topic  in  itself  (3,  1J),  14)  and 
is  currently  considered  the  weak  link  in 
error  analysis  of  the  factor  of  safety 
value.  Spatial  variability  could  be  as- 
sessed through  three  components:  (1)  the 
natural  heterogeneity  of  the  soil, 
(2)  measurement  error,  and  (3)  limited 
information  about  subsurface  conditions. 
Anderson  (2^)  discusses  the  merits  of  con- 
sidering model  error  when  assessing  the 
variability  of  the  factor  of  safety. 


DISTRIBUTION  PROPERTIES 


Data  were  collected  at  two  sites — 
a  silver-lead-zinc  mine  in  the  Coeur 
d'Alene  mining  district  of  Idaho  (CDA 
data),  and  a  Southwest  U.S.  copper  mine 
(SW  data).  The  appendix  describes  the 
field  collection  of  the  data.  Table  1 
enumerates  basic  statistics  from  these 
two  sites;  effective  stress  is  used 
throughout. 

Note  that  a  rather  large  number  of  sam- 
ples were  obtained  at  each  site,  giving  a 


significant  population  size  for  consider- 
ation. The  variability  in  the  internal 
angle  of  friction  (<))),  cohesion  (c),  and 
soil  density  (y)  is  consistent  with  the 
results  of  Lee  and  Singh's  (11)  earlier 
works.  In  fact,  <(>  has  less  variation 
than  anticipated,  although  Singh's  work 
focused  on  variability  between  laborato- 
ries rather  than  on  within-laboratory 
variability.  Note  the  close  relationship 
between  the  values  of  Bishop's  factor  of 
safety  and  that   of  Fellenius.   Since  the 


propagation  of  error  work  uses  Fellenius' 
factor  of  safety,  we  will  only  use  Felle- 
nius' numbers  henceforth. 

To  better  visualize  the  data,  histo- 
grams of  the  data  are  shown  in  figures  3 
through  6.  In  figure  3,  histogram  of  the 
factor  of  safety,  note  the  wide  variabil- 
ity in  the  CDA  data.  Over  10  pet  of  the 
data  lies  to  the  left  of  1;  thus  the 
probability  of  failure  is  predicted  to 
be  quite  high.  However,  there  are  two 
points  to  keep  in  mind: 


paper  considers  the  normal,  log-normal, 
and  Weibull  distributions.  Table  2  sum- 
marizes curve-fitting  these  functions  to 
Fellenius'  factor  of  safety  values. 


Two  interesting  points 
from  table  2: 


may  be  drawn 


1 .  The  logarithmic  transformation  does 
not  improve  fitting  the  PDF's.  This  is 
reasonable,  for  the  log-normal  transfor- 
mation works  best  for  data  that  are  ini- 
tially skewed  to  the  left. 


1.  In  this  case,  the  material  is  con- 
tained behind  a  berm  of  borrow  material; 
thus  the  factor  of  safety  values  are  sim- 
ulated, not  actual  values. 

2.  This  wide  variability  gives  credi- 
bility to  the  variance  reduction  concept, 
which  is  intrinsic  to  the  probabilistic 
work  of  Vanmarcke  and  Anderson  (2^,  16) . 
Thus,  variance  reduction,  although  a  con- 
troversial concept  as  a  "pragmatic  value" 
at  this  time,  is  conceptually  acceptable 
and  at  least  partially  verified  as  to 
need  by  this  distribution. 

Certainly  the  histogram  for  the  factor 
of  safety  (fig.  3)  demonstrates  a  strong 
central  tendency  in  both  cases;  the  next 
logical  question  is  whether  a  conven- 
tional mathematical  PDF  "fits"  the  data 
properly  or  not.  There  are  many  candi- 
date functions  for  consideration;   this 


2.  The  Weibull  PDF  does  not  improve  on 
the  curve  fits  for  the  normal  PDF.  Since 
the  Weibull  is  more  "obscure"  (for  exam- 
ple, parameters  are  not  easily  under- 
stood, and  the  curve-fitting  is  appreci- 
ably more  difficult) ,  there  seems  to  be 
no  reason  to  use  it  in  preference  to  the 
normal  PDF. 

It  is  true  that  there  are  other  candi- 
date distributions — the  beta,  gamma, 
Poisson,  et  al. — that  could  have  been 
considered.  In  particular,  the  geotech- 
nical  community  has  shown  some  preference 
for  the  beta  PDF.  Upon  investigation  of 
the  properties  of  the  Weibull  and  beta 
PDF's,  it  is  found  that  both  are  versa- 
tile functions  that  can  fit  a  variety  of 
situations  in  a  very  similar  manner.  Our 
choice  is  to  investigate  just  the  Weibull 
PDF  in  this  report.  It  is  felt  that  this 
investigation  is   sufficient   and  the 


TABLE  2.  -  Three  model  curve  fits  to  Fellenius'  factor  of  safety 


Probability 

density 
functions 1 

Index  of 
determination2 

X2  values3 

Probability 

density 
functions 1 

Index  of 
determination2 

X2  values3 

CDA  SITE 

SW  SITE 

Log  normal. . . 

0.62 
.59 
.62 

418.41 
NA 
NA 

Log  normal . . . 

0.96 
.94 
.96 

56.13 

NA 
NA 

NA  Not  available. 

1A11  curve  fitting  done  using  Gauss'  nonlinear  iterative  technique  (13). 


Measures  how  well  the  data  agree  with  curve-fit  function.   Reduces  to  the  correla- 

used  for  comparison  is  9.49  for  these 


tion  coefficient  in  the  linear  case  (13) . 
3Goodness-of-fit  test;   the  critical  value 


cases. 

4 


Reject. 
5 Accept. 


20 


>         0 

o 

Z       40 

M 

O 


3d 


Coeur    d'    Alene 


20 


■  : 


Southwest 


F  =  1.95 
Sc=0.24 


0.6  1.0  1.4  1.8  2.2  2.6  3.0 

F    (FELLENIUS) 

FIGURE  3.  -   Factor  of  safety  histograms. 


I 

I      I 

^^r\         Coeur 

flT 

Alene 

10 

c   =3.53 
Sc=2.94 



0 

H 

i    Vr 

I 

40 

> 
O 

z 

M 

O      30 

LU 
X 


20 


10 


COHESION,    psi 

FIGURE  5.  -  Cohesion   histograms. 


30 


20 


10 


40 


30 


20 


10 


25       30       35       40       45       50 
PHI  ANGLE,  deg 

FIGURE  4.  -   Internal  friction  angle  histograms. 


20 


uj  0 

°      30 


Coeur    d'    Alene 


Southwest 


FIGURE  6.  -  Soil  density  histograms. 


normal  PDF  is  as  good  a  choice  as  can  be 
made  for  these  two  examples.  It  is  true 
that  the  normal  distribution's  curve  fit 
to  the  CDA  is  not  outstanding  as  indi- 
cated by  the  chi-square  measure.  The 
data  are  simply  not  conducive  to  a  good 
curve  fit,  regardless  of  what  function 
one  might  try. 

Table  3  enumerates  pertinent  statistics 
and  curve-fit  parameters  for  internal 
angle  of  friction,  cohesion,  and  density. 

In  view  of  the  preceding  discussion, 
and  since  the  remaining  PDF's  are  of 
lesser  interest  to  us,  only  the  normal 
curve  will  be  used  to  fit  these  histo- 
grams. The  internal  angle  of  friction 
distributions  (fig.  4),  although  again 
demonstrating  central  tendency,  are  each 
skewed,  both  to  the  right.  The  index  of 
determination  for  both  fits  is  satisfac- 
tory, with  a  somewhat  better  fit  for  the 
SW  data.  Note  that  the  fitted  parameters 
vary  somewhat  from  the  sample  means  and 
standard  deviations;  the  degree  of  vari- 
ability is  a  function  of  the  goodness  of 
fit. 


It  is  not  surprising  that  the  cohesion 
value  (fig.  5)  has  the  largest  variabil- 
ity, for  this  is  apparently  inherent  in 
the  cohesion  parameter  (11) .  Thus  the 
normal  curve  approximations  to  the  cohe- 
sion histograms  are  not  particularly  im- 
pressive, as  indicated  by  the  relatively 
low  indices  of  determination.  Here  the 
histograms  would  indicate  a  slight  skew- 
ness  to  the  left. 

Finally,  consider  the  histogram  (fig. 
6)  for  the  soil  density.  Again,  curve 
fitting  these  curves  is  of  somewhat  nebu- 
lous value;  for  example,  a  uniform  dis- 
tribution may  be  a  better  candidate  for 
the  SW  data.  But  to  be  consistent,  the 
normal  curve  approximation  was  used;  the 
resulting  indices  of  dispersion  reflect 
the  degree  of  disagreement  between  the 
data  and  fitted  function. 

To  summarize,  central  tendency  is  ex- 
hibited by  all  of  the  histograms  (and 
corresponding  PDF's)  considered.  The 
normal  curve  seems  to  be  an  acceptable 
candidate  model.  For  consistency's  sake, 
the  normal  curve  was  then  used  for  all 


TABLE  3.  -  Normal  curve-fit  parameters  for  <j> ,  c, 
y,  and  F 


Normal  curve-fit 


parameter 


Area 


Mean 


Index  of 
determination2 


CDA  SITE 


♦ 

c 

Y 

F-Fellenius 


.deg.. 
.psi. . 
.pcf. . 


118.78 

115.22 

147.38 

14.14 


0.82 
.79 
.62 
.62 


SW  SITE 


* : 

c 

Y 

F-Fellenius, 


,deg.. 
.psi. . 
.pcf. . 


291.09 

296.44 

441.30 

17.72 


.91 
.85 
.88 
.96 


1  Equation  for  normal  PDF  used: 

-  (x  -  mean)  2/2  (std.  dev.)2 


y  = 


(area) 


/2rr  (std.  dev.) 
where   y  =  dependent  variable  (histogram  value) , 
x  =  independent  variable  (F,  <)) ,  c,  y)> 
2A  measure  of  goodness  of  fit  for  nonlinear  functions; 
see  reference  13  for  details. 


histogram  curve  fitting.  Also,  the  se- 
verity of  concern  decreases  considerably 
as  one  moves  from  the  factor  of  safe- 
ty to  the  supporting  data  sets.  Thus, 
for  probabilistic  work,   in  finding  the 


probability  of  failure  and  in  setting 
confidence  intervals  on  the  estimated 
mean  factor  of  safety  value,  these  two 
examples  fortify  use  of  the  normal  proba- 
bility density  function. 


PROPAGATION  OF  ERROR  APPLIED  TO  THE  DATA 


As  stated  previously,  it  is  desirable 
to  take  only  a  few  field  samples  and  use 
the  propagation  of  error  technique  (6)  to 
estimate  the  variance  of  the  factor  of 
safety.  However,  since  we  have  complete 
sets  of  data  available,  we  naturally  used 
all  these  data  for  comparative  slope  sta- 
bility work  herein.  (See  appendix  for 
the  data  description.)  Table  4  shows  the 
coefficient  of  variation  for  Fellenius' 
factor  of  safety  value  with  assumed  pore 
pressure  coefficient  of  variation  values 
of  0,  10,  and  20  pet.  These  values  are 
compared  to  the  sample  coefficient  of 
variation  value.  In  the  case  of  the  SW 
embankment ,  the  agreement  between  the 
propagation  values  and  that  of  the  sample 
data  is  quite  good.  The  CDA  data  do  not 
agree  as  well,  but  the  divergence  would 
seem  acceptable.  Note  that  the  pairwise 
correlations  between  the  soil  parameters 
(for  example,  "<j>"  and  "c")  may  or  may  not 
be  included  in  the  computations,  depend- 
ing on  the  user's  preference.  We  chose 
to  use  these  correlation  coefficients  in 
our  computations.   Table  5  enumerates  the 


pairwise  correlation  coefficients  used  in 
calculating  table  4. 

TABLE  4.  -  Comparison  of  coefficient 
of  variation  estimates  for  Felle- 
nius' factor  of  safety 


Coefficient  of  variation 


0  pet1 

10  pet1 

20  pet1 

Actual  sample  values2 


1The  entries  are  the  propagation  of  er- 
ror computed  values  for  the  specified 
pure  pressure. 

2 The  entries  are  the  comparative  values 
obtained  from  the  available  sample  data — 
see  table  1. 

TABLE  5.  -  Correlation  coefficients 


Correlation  coefficient 
rc,4> 

rc,Y 

r<J>,Y 


CONFIDENCE  INTERVAL  FOR  THE  FACTOR  OF  SAFETY 


There  is  one  basic  statistical  approach 
that  can  be  easily  applied  to  the  factor 
of  safety  computation  which  provides  a 
simplistic  marriage  between  the  world  of 
probability  and  the  presently  used  deter- 
ministic computation — putting  a  confi- 
dence limit  about  the  factor  of  safe- 
ty value.  A  brief  review  of  statistical 
principles  would  seem  in  order.  Suppose 
one  is  concerned  with  a  variable  whose 
PDF  is  not  necessarily  normally  distri- 
buted, but  for  which  the  population  mean 
and  standard  deviation  exist.  Next,  a 
sample  size  "n"  is  chosen.  By  some  ran- 
dom process,  "n"  samples  are  selected 
from  the  original  population  and  a  sample 
mean  (x)  is  computed.  This  process  is 
then  repeated  until   a   population  of 


sample  means  (all  of  which  are  of  the 
preselected  sample  size  "n")  exists.  The 
Central  Limit  Theory  of  Statistics  states 
that  the  PDF  of  mean  values  will  be  nor- 
mally distributed  provided  "n"  is  suffi- 
ciently large.  Furthermore,  the  popula- 
tion mean  of  the  new  distribution — the 
distribution  of  sample  means — is  the  same 
as  the  population  mean  of  the  original 
PDF,  but  the  standard  deviation  of  the 
new  distribution  equals  the  original 
PDF's  standard  deviation  divided  by  the 
square  root  of  "n."  When  the  Central 
Limit  Theorem  is  actually  applied  to  a 
data  set,  only  one  sample  of  size  "n"  is 
taken  from  which  the  sample  mean  and 
standard  deviation  are  computed. 
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This  powerful  theorem  can  be  applied  to 
the  distribution  of  factor  of  safety  val- 
ues. The  key  of  the  theorem  is,  of 
course,  how  many  samples  are  required  for 
it  to  be  valid.  The  answer  to  this  query 
is  certainly  not  simple;  if  the  original 
distribution  is  very  similar  to  the  nor- 
mal curve,  then  small  sample  sizes  suf- 
fice. For  symmetrical  distributions  (the 
preceding  work  and  that  of  Baecher  and 
Marr  (_3)  substantiate  that  one  may  assume 
symmetry) ,  a  sample  size  of  eight  or 
greater  would  seem  appropriate.  Thus  the 
computational  procedure  would  be  as 
follows: 


1.  Collect  a  representative  set  of 
soil  properties;  a  minimum  of  eight  sets 
is  recommended.  Compute  a  corresponding 
set  of  factor  of  safety  values. 

2.  Compute  the  sample  mean  (F)  and 
standard  deviation  (sF)  for  these  factor 
of  safety  values. 

3.  Apply  the  Central  Limit  Theorem  by 
forming  a  95-pct  confidence  interval  for 
the  true  but  unknown  population  mean  fac- 
tor of  safety  value  (up).  The  following 
computation  is  requisite: 


where 


F  ±  to. 975  (n~l)  *  sf/  ^i  =  95-pct  confidence  interval, 
n  =  sample  size, 


and   tg. 975  (n-1)   =  "t"  distribution  value,  which  depends  on  sample  size 

and  the  level  of  confidence. 


-  1 


There  is  no  set  rule  that  mandates  a  95- 
pct  confidence  interval;  this  is  a  con- 
ventional and  conservative  choice. 

Let  us  investigate  what  it  meant  by  a 
confidence  interval.  For  example,  a  95- 
pct  confidence  interval  implies  that  if 
one  were  to  form  100  such  intervals , 
95  of  these  would  contain  the  true  (but 
unknown)  population  parameter.  In  prac- 
tice, one  computes  only  one  such  interval 
and  relies  on  statistical  theory  to  pro- 
vide the  statement  of  confidence  or 
assurance. 

A  simulation  may  prove  illustrative. 
The  authors  selected,  via  random  sam- 
pling, 16  sample  factor  of  safety  values 
from  each  of  the  previously  discussed 
data  populations.  It  was  assumed  that 
the  computed  sample  mean  value  was  a  best 
approximation  to  the  true  factor  of  safe- 
ty value.  Then  confidence  intervals  for 
subsamples  of  size  4  were  formed,  then  of 
size  8,   and  finally  of  size  16.   These 


confidence  intervals  are  shown  in  figures 
7  and  8  for  the  two  data  populations  of 
interest.  Note  that  in  all  but  one  case 
the  95-pct  confidence  intervals  did  con- 
tain the  overall  sample  mean.  Also,  note 
that  the  intervals  tend  to  lessen  in 
length  as  the  sample  size  increases.  In 
figure  7 — the  CDA  data — samples  of  size  4 
include  factor  of  safety  values  less  than 
1  in  three  of  the  four  cases.  Also,  one 
of  the  two  samples  of  size  8  extends  be- 
low a  factor  of  safety  value  of  1.  This 
is,  of  course,  an  indicator  of  concern 
for  the  embankment  designer. 

Again,  the  preceding  simulation  is  tu- 
torial, for  in  practice  one  would  obtain 
only  one  confidence  interval  and  use  it 
to  make  a  statement  about  the  factor  of 
safety.  Furthermore,  an  example  is,  at 
best,  illustrative  and  certainly  does 
not  replace  theoretical  considerations. 
Ideally,  the  interval  would  contain  no 
factor  of  safety  values  below  the  criti- 
cal value  of  1. 


GEOTECHNICAL  INNOVATION  AND  PROBABILISTIC  MODELING 


There  is  a  new  approach  regarding 
the  "best"  model  for  the  factor  of 
safety   calculation.   The   traditional 


two-dimensional  slip  circle  of  Fellenius 
and  Bishop  is  being  challenged.  In  par- 
ticular, a  three-dimensional  cylinder  has 
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FIGURE  7.  -  Ninety-five  percent  confidence 
interval-factor  of  safety,  CDA  Mine— with  vary- 
ing sample  size. 

been  postulated  by  Vanmarcke  and  others 
(—>  16) •  Theoretical  development  reduces 
the  cylindrical  model  effectively  to  that 
of  Bishop  (5,  9_) ,  but  with  end-effect 
forces  considered. 

A  most  important  concept  in  the  geo- 
technical  works  of  Vanmarcke  (16)  and 
Anderson  (2^)  is  variance  reduction.  This 
rather  complex  topic  is  explained  in 
depth  in  the  two  aforementioned  refer- 
ences; the  following  is  a  heuristic  over- 
view. Rather  than  being  interested  in 
the  usual  statistical  value  of  the  vari- 
ance of  the  strength  parameter,  the  geo- 
technical  engineer  is  concerned  when  a 
succession  of  weak  strength  values  oc- 
curs. In  such  a  case,  the  possibility 
of  failure  is  greatest.  If  one  were  to 
envision  a  moving  average  curve  being  fit 
to  a  spatial  succession  of  strength  val- 
ues,  the   variability   from   the   moving 
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FIGURE  8.  -  Ninety-five  percent  confidence 
interval— factor  of  safety,  SW  Mine— with  varying 
sample  size. 

averages  is  representative  of  the  desired 
reduced  variance  (fig.  9).  A  scaling 
factor,  variance  reduction,  is  used  to 
reduce  the  statistical  variance  value  to 
the  desired  variance  value.  The  actual 
computation  of  the  scaling  factor  is  not 
as  simple  as  described  herein.  Involved 
is  a  plot  of  the  autocorrection  function 
for  the  values  of  concern.  Next,  a  curve 
is  fit  to  the  autocorrelation  function, 
and  said  function  is  included  in  an  inte- 
gration that  computes  the  desired  reduced 
variance  value. 

Two  major  problems  in  the  application 
of  the  probabilistic  model  were  found 
during  a  recent  field  test  under  a  Bureau 
of  Mines  contract.  One  concerns  the  ef- 
fects of  variance  reduction  on  the  prob- 
ability of  failure,  and  the  other  is  the 
high  cost  of  obtaining  strength  data 
using  a  penetrometer. 

The  variance  reduction  value  proved  to 
be  quite  severe,  effectively  reducing  the 
factor  of  safety  PDF  to  a  single  spike. 
The  dilemma  is  that  one  has  gone  through 
the  complex  mechanism  of  probabilistic 
slope  stability  and  has  only  obtained  a 
univalued  answer,  as  in  the  deterministic 
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FIGURE  9. 


DISTANCE    ALONG    EMBANKMENT    CENTERLINE 

Variance  reduction  applied  to  spatially  oriented  shear  strength  values. 


approach.  Only  replicated  fieldwork  will 
verify  whether  or  not  this  is  the  general 
case. 

To  accommodate  the  computations  for 
variance  reduction,  a  sequence  of  closely 
spaced  strength  values  must  be  obtained 
in  all  three  dimensions  of   the  dam  or 


embankment.  Cone  penetrometer  soundings 
are  required  to  gather  the  necessary  in- 
formation. This  adds  greatly  to  the  cost 
of  obtaining  the  factor  of  safety.  The 
obvious  question  to  be  answered  is  wheth- 
er the  benefits  gained  from  a  probabilis- 
tic analysis  justify  the  increased  cost. 


CONCLUSIONS 


A  large  sample  of  direct  shear  tests 
was  taken  at  two  waste  embankment  sites. 
Of  the  statistical  model  tested,  the  nor- 
mal probability  density  function  best  fit 
these  data.  This  conclusion  is  site  spe- 
cific and  should  not  be  imprudently  ex- 
trapolated. The  authors  are  aware  that  a 
wealth  of  data  is  now  available  to  the 
mining  industry,  since  penetrometers  of 
various  kinds  are  now  being  used  more 
frequently.  Thus,  detailed  investigation 
of  the  variation  of  strength  with  depth 
can  now  be  routinely  incorporated  to  make 
computations  more  meaningful.  In  fact, 
the  Bureau  has  contracted  for  soil  param- 
eter collecting  using  cone  data  ( 3) .  The 
Bishop  code  with  the  propagation  of  error 


and  probability  of  failure  calculations 
is  available  from  the  Bureau  for  inter- 
ested users. 

A  suggested  first  step  towards  intro- 
ducing useful  probabilistic  concepts  into 
the  factor  of  safety  computation  is  to 
compute  a  confidence  interval  about  the 
average  factor  of  safety  value ,  replacing 
current  usage  of  a  conservative  determin- 
istic factor  of  safety  value.  It  would 
seem  most  advantageous  not  only  to  know 
the  factor  of  safety  value  but  to  have  a 
sense  of  the  error  in  that  value  and  an 
indication  of  how  good  or  reliable  the 
factor  of  safety  might  be. 
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APPENDIX. —DATA  COLLECTION,  SLOPE  STABILITY,  AND  LABORATORY  TESTING 


DATA  COLLECTION--CDA 

A  detailed  investigation  of  the  im- 
poundment facility  was  not  undertaken, 
since  the  scope  of  this  project  was  to 
determine  strength  variability  in  the 
beach  material  and  to  model  an  embankment 
consisting  solely  of  this  material. 

The  embankment  was  constructed  by  the 
upstream  method  using  a  starter  dike  of 
borrow  material.  The  dike  is  built 
against  a  steep  side  of  a  valley  with  an 
overall  downstream  slope  of  2:1.  The  im- 
poundment area  covers  approximately  7 
acres.  Tailings  were  deposited  from 
spigot  lines  extended  from  the  crest  of 
the  embankment.  They  contain  from  15  to 
50  pet  material  passing  the  No.  200  U.S. 
Standard  Sieve.  In  general,  the  coarser 
material  settled  out  close  to  the  crest, 
but  a  layering  effect  became  evident 
while  sampling.  Fluctuating  pond  size 
and  different  spigot  locations  could  have 
caused  this  phenomenon. 


for  the  computer  model.  These  were  se- 
lected to  evaluate  a  potentially  unsafe 
condition  (that  is,  P  (failure)  >  0)  and 
in  no  way  reflect  the  actual  conditions 
of  the  embankment  from  which  the  soil  was 
extracted.  The  location  of  the  phreatic 
surface,  based  on  the  10  pet  freeboard 
headwater,  was  determined  by  the  finite- 
element  method  (4).2  No  tail  water  was 
assumed  in  the  model.  The  material  in 
the  embankment  was  assumed  to  be  homoge- 
neous to  the  extent  that  no  layers  were 
coded  into  the  model.  The  code  was  exe- 
cuted using  74  sets  of  strength  parame- 
ters (c,  <j) ,  y,)  in  order  to  calculate  the 
coefficent  of  variation  of  the  factor  of 
safety  and  compare  it  to  the  value  ob- 
tained by  using  the  propagation  of  error 
formula.  The  value  for  density  that  was 
used  in  each  run  was  an  average  of  the 
densities  of  the  samples  used  to  conduct 
the  direct  shear  test  for  a  corresponding 
c,  <J>. 

DATA  COLLECTION— SW 


Data  were  taken  along  the  perimeter  of 
the  embankment  15  ft  from  the  crest  at 
intervals  of  150  ft.  At  each  data  loca- 
tion, four  thin-wall,  3-in-diam  Shelby1 
tube  samples  were  taken  at  the  corners  of 
an  imaginary  square  measuring  2  ft  on 
each  side.  Before  inserting  the  Shelby 
tubes  by  hand,  the  tailings  were  in- 
spected for  disturbance  and  cracking. 
If  either  of  these  conditions  existed, 
the  cluster  of  four  tubes  was  displaced 
slightly  to  avoid  the  irregularity  or 
cracks.  When  the  tubes  were  extracted, 
the  ends  were  capped  and  sealed,  and  the 
tubes  were  placed  in  a  foam-lined  con- 
tainer for  transportation  by  truck  to  the 
laboratory.  Care  was  taken  that  minimum 
sample  disturbance  occurred. 

STABILITY  ANALYSIS— CDA 

A  slope  of  2  to  1  and  a  phreatic  sur- 
face having  10  pet  freeboard  were  chosen 

1 Reference  to  specific  trade  names  is 
made  for  identification  only  and  does  not 
imply  endorsement  by  the  Bureau  of  Mines. 


As  with  the  CDA  embankment,  an  exten- 
sive field  investigation  was  not  under- 
taken. A  brief,  general  description  will 
again  be  presented. 

The  SW  tailings  embankment  has  been  in- 
active for  years.  It  was  constructed 
from  copper  tailings ,  using  the  upstream 
spigoting  method.  The  slurry  was  about 
45  pet  tailings  by  weight;  the  tailings 
grind  was  40  pet  passing  the  No.  200 
Sieve  and  75  pet  passing  the  No.  65 
Sieve.  The  embankment  is  bounded  on  one 
side  by  another  tailings  embankment.  The 
tailings  pond  covers  approximately  300 
acres  and  averages  300  ft  high.  The 
overall  slope  is  approximately  2  to  1. 

The  data  for  the  SW  embankment  were 
also  collected  using  thin-wall  Shelby 
tubes.  The  embankment  had  from  6  in  to 
several  feet  of  overburden  placed  on 
top  of   the  tailings   for  dust   control 

^Underlined  numbers  in  parentheses  re- 
fer to  items  in  the  list  of  references 
preceding  the  appendix. 


15 


purposes.  The  top  3  to  5  ft  of  tailings 
was  compacted  from  equipment  placing  the 
overburden;  in  some  cases,  desiccation 
was  apparent.  An  attempt  was  made  to  in- 
sert the  Shelby  tubes  using  a  jack,  with 
resistance  provided  by  a  truck.  This 
caused  too  much  sample  disturbance  so  a 
front-end  loader  was  employed  to  scrape 
off  the  overburden  and  push  the  Shelby 
tubes  into  the  embankment.  (The  tubes 
inserted  by  the  jacks  were  not  included 
in  the  data  analysis.) 


SAMPLE  TESTING 

Standard  direct  shear  tests  were  per- 
formed on  a  series  of  samples  extracted 
from  the  Shelby  tubes  with  a  circular 
trimmer  ring.  Owing  to  the  time  needed 
to  test  all  the  soil  samples,  some  corro- 
sion took  place  inside  the  Shelby  tubes. 
These  samples  were  sawed  in  half  to  re- 
duce sample  disturbance  when  the  soil  was 
extracted  from  the  tube  into  the  trimmer 
ring. 


The  tubes  were  extracted  by  hand, 
sealed,  and  packed  into  foam-lined  boxes 
for  minimal  disturbance.  It  was  neces- 
sary to  lay  the  tubes  in  a  horizontal  po- 
sition for  shipping. 

The  data  station  pattern  is  shown  in 
figure  A-l.  This  pattern  was  chosen  so 
that  correlations  might  be  studied  based 
on  distance  from  the  crest  and  spacing 
distance  between  data  points. 

Unlike  the  CDA  embankment,  the  computer 
model  was  set  up  to  reflect  the  geometry 
and  phreatic  surface  location  of  the  ac- 
tual embankment.  The  cross  section  that 
was  used  is  shown  in  figure  A-2.  The  ac- 
tual embankment  had  a  starter  dike  that 
was  not  included  in  the  analysis.  This 
exclusion  could  produce  lower  factors  of 
safety  than  actually  exist.  All  130 
triplets  (c,  <f> ,  y,)  were  run  through  the 
computer  code  to  obtain  the  130  values 
of  the  factor  of  safety  necessary  to 
calculate  the  PDF  and  coefficient  of 
variation. 


Each  sample  was  consolidated  before 
shearing  by  applying  normal  load  to  allow 
for  drainage.  Initial  displacement  mea- 
surements from  consolidation,  along  with 
direct  shear  test  data,  are  given  in  ta- 
bles A-l  and  A-2.  At  a  minimum,  normal 
stresses  of  25  psi,  50  psi,  and  100  psi 
were  used.  For  the  CDA  data,  an  addi- 
tional point  of  75  psi  was  used  to  deter- 
mine each  <f> ,  c  combination.  Owing  to  the 
layering  effects  of  each  embankment ,  it 
was  necessary  to  repeat  various  normal 
loads  to  ensure  that  enough  data  were 
available  for  a  good  least  squares  curve 
fit  for  <|>  and  c.  Average  shearing  rates 
for  the  tests  were  0.305  mm/sec  for 
the  25-psi  normal  load,  0.275  mm/sec  for 
50  psi  and  75  psi,  and  0.25  mm/sec  for 
100  psi.  The  slow  rate  was  used  to  pre- 
vent pore  pressure  buildup  and  ensure 
drainage. 

Cumulative  grain-size  distribution 
tests  were  run  for  14  tubes  from  the  CDA 
data.  At  least  one  tube  was  tested  from 
each  data  location  in  an  attempt   to 


•   • 


•  •    •   •    • 

•  •    •    •    • 
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324' 


Embankment!   crest 


FIGURE  A-l.  -   SW  data  station  pattern 


16 


represent  the  range  of  data  collected.  A 
range  of  distributions  is  shown  in  figure 
A-3. 

COMPUTER  PROGRAM  FOR  PROPAGATION 
OF  ERROR 


The  computer  program  used  to 
the  coefficient  of  variation  of 
tor  of  safety  is  a  modified  vers 
program  written  by  Bailey  (4). 
gram  uses  Bishop's  Simplified 
Slices  to  determine  the  minimum 
safety  for  a  given  embankment 
The  failure  surface  is  assumed 
arc  of  a  circle. 


calculate 
the  fac- 
ion  of  a 
The  pro- 
Method  of 
factor  of 
geometry, 
to  be  an 


The  program  divides  the  slope  cross 
section  into  vertical  slices  and  evalu- 
ates moments  about  the  circle  center 
using  the  force  vectors  at  the  base  of 
each  slice.  The  horizontal  side  forces 
and  vertical  shear  forces  on  each  slice 
are  not  used  in  the  calculations.  The 
number  of  slices  is  input  by  the  user 
but  may  be  altered  slightly  by  the 
program  so  that  slice  boundaries  coin- 
cide with  selected  points  on  the 
embankment . 


500 


The  slice  boundaries  are  determined  be- 
fore any  trial  circles  are  analyzed,  so 
the  number  of  slices  in  each  sliding  mass 
varies  with  the  radius  of  the  trial  cir- 
cle.  To  determine  a  minimum  factor  of 
safety,  the  program  evaluates  1,331  trial 
circles.   The  center  of  each  circle  is  a 
point  on  an  imaginary  square  grid  placed 
above  the  slope  of  the  embankment.   At 
each  point  on  the  control  grid,   the  pro- 
gram computes  a  minimum  radius;   that  is, 
a  radius  that  will  just  touch  the  slope 
and  a  maximum  radius  based  on  the  geome- 
try of  the  slope  profile.   A  series  of 
factors  of  safety  is  computed,  starting 
with  a  radius  slightly  smaller  than  the 
maximum  radius.   Each  successive  radius 
is  reduced  by  approximately  one-eleventh 
of  the  difference  between  the  minimum  and 
maximum  radii. 

Coefficients  for  evaluation  of  equation 
4  are  required  for  input  in  addition  to 
embankment  geometry,  phreatic  surface  in- 
formation, and  soil  parameters.  The  pro- 
gram outputs  the  coefficient  of  vari- 
ation and  standard  deviation  for  the  min- 
imum factor  of  safety,  the  probability  of 
failure  (assuming  normal  statistics),  and 


200 
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HORIZONTAL  DISTANCE,   ft 
FIGURE  A-2.  -   Idealized  profile-SW  embankment. 
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the  contribution  by  term  (equation  4)  to 
the  standard  deviation  of  the  factor  of 
safety. 

The  pore  pressure  in  the  computer  pro- 
gram is  calculated  at  the  base  of  each 
slice,  based  on  static  head.  Similarly, 
the  contribution  to  the  error  term  by 
pore  pressure  is  on  a  slice-by-slice 
basis.  Since  the  user  inputs  the  coef- 
ficients  of  variation,  it   is  convenient 


to  calculate  both  the  variance  term 
and  the  covariance  terms  for  pore  pres- 
sure. In  our  calculations,  it  was  as- 
sumed that  the  correlations  between  pore 
pressures  and  the  strength  parameters 
=  0,  and  that  the  correlations  between 
pore  pressures  for  any  two  slices  =  1. 
The  covariance  terms  for  the  strength 
parameters   are  included;  that  is,   rc  x 

rc,Y,  and  r<J>,Y. 
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TABLE  A-2.  -  SW  Laboratory  data 


Fellenius 

Wet 

Internal 

Cohesion, 

Fellenius 

Wet 

Internal 

Cohesion, 

Sample 

factor  of 

density, 

friction 

psi 

Sample 

factor  of 

density, 

friction 

psi 

safety 

pcf 

angle,  deg 

safety 

pcf 

angle,  deg 

1 

1.948 

111.36 

39.39 

4.22 

66 

2.107 

100.08 

43.20 

2.60 

2 

1.932 

103.54 

39.44 

3.96 

67 

1.875 

114.14 

38.10 

2.30 

3 

2.124 

103.67 

38.51 

9.16 

68 

1.674 

105.11 

39.50 

.00 

4 

1.899 

115.81 

37.13 

2.40 

69 

1.904 

108.95 

36.31 

7.34 

5 

1.838 

112.50 

36.25 

2.40 

70 

1.775 

111.33 

40.98 

.05 

6 

2.008 

105.31 

36.80 

5.52 

71 

1.871 

108.26 

39.54 

2.40 

7 

1.969 

95.74 

36.42 

5.26 

72 

1.801 

122.41 

36.60 

2.60 

8 

2.137 

100.14 

44.04 

.57 

73 

1.790 

119.02 

38.60 

1.20 

9 

2.017 

103.96 

37.51 

4.74 

74 

1.964 

107.91 

41.00 

2.40 

10 

2.245 

99.51 

41.65 

3.70 

75 

1.819 

107.18 

34.80 

7.30 

11 

2.402 

102.52 

42.03 

6.30 

76 

1.913 

110.10 

39.80 

1.60 

12 

1.937 

111.54 

43.19 

.00 

77 

1.913 

105.47 

40.40 

1.20 

13 

2.009 

96.73 

36.20 

6.30 

78 

2.040 

112.15 

41.90 

2.60 

14 

1.713 

105.49 

31.90 

5.50 

79 

1.887 

110.95 

40.80 

.60 

15 

1.550 

104.15 

36.91 

.00 

80 

1.661 

114.64 

35.20 

3.30 

16 

2.246 

104.24 

43.23 

1.61 

81 

1.622 

109.11 

34.90 

2.90 

17 

2.388 

106.46 

42.78 

4.74 

82 

2.363 

119.12 

43.80 

6.90 

18 

1.820 

103.97 

34.56 

4.48 

83 

2.288 

112.54 

47.00 

.50 

19 

2.328 

109.54 

37.78 

11.25 

84 

2.030 

114.79 

45.00 

.00 

20 

2.609 

103.14 

43.77 

13.59 

85 

2.102 

106.93 

42.40 

3.40 

21 

1.909 

110.85 

39.89 

2.66 

86 

1.916 

103.65 

40.14 

1.35 

22 

2.055 

104.04 

42.71 

2.05 

87 

1.984 

104.89 

40.24 

3.96 

23 

2.169 

97.81 

44.12 

2.66 

88 

1.962 

102.07 

39.95 

3.96 

24 

2.216 

103.96 

40.59 

8.38 

89 

2.524 

113.11 

43.23 

11.51 

25 

2.134 

100.61 

42.96 

1.61 

90 

2.144 

109.08 

37.78 

10.73 

26 

2.075 

102.06 

42.82 

2.40 

91 

2.096 

112.27 

36.69 

11.25 

27 

1.900 

108.50 

43.10 

.00 

92 

1.771 

98.86 

41.10 

.00 

28 

1.896 

99.67 

36.28 

7.99 

93 

1.980 

100.70 

41.90 

1.80 

29 

1.812 

96.12 

39.59 

1.61 

94 

2.041 

100.83 

44.10 

.30 

30 

1.965 

105.12 

41.27 

2.14 

95 

1.802 

105.89 

41.60 

.00 

31 

2.617 

111.45 

42.90 

15.50 

96 

1.987 

107.55 

40.80 

3.20 

32 

2.194 

112.55 

44.50 

.90 

97 

2.126 

111.88 

35.60 

14.80 

33 

1.802 

104.27 

34.06 

7.73 

98 

1.940 

111.79 

43.70 

.00 

34 

1.998 

115.73 

37.82 

7.86 

99 

1.928 

100.15 

40.90 

2.10 

35 

1.981 

108.40 

35.98 

9.42 

100 

1.974 

111.59 

37.50 

7.40 

36 

1.773 

103.94 

41.13 

.00 

101 

2.212 

102.85 

43.86 

4.70 

37 

2.044 

118.10 

33.80 

15.67 

102 

2.204 

97.90 

42.91 

5.00 

38 

2.150 

112.20 

42.27 

2.40 

103 

2.110 

99.18 

41.89 

4.48 

39 

1.740 

109.56 

33.15 

7.60 

104 

1.894 

105.49 

37.56 

5.52 

40 

1.716 

114.06 

40.20 

.00 

105 

1.945 

98.80 

41.17 

2.14 

41 

2.360 

104.86 

40.60 

12.80 

106 

1.836 

111.97 

27.10 

18.79 

42 

1.602 

109.29 

31.00 

7.00 

107 

1.923 

102.53 

31.00 

15.30 

43 

1.971 

115.27 

41.30 

1.30 

108 

1.574 

111.19 

31.80 

5.50 

44 

2.028 

101.66 

39.80 

5.50 

109 

1.960 

101.19 

36.10 
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